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Abstract. Astrophysical turbulence is magnetohydrodynamic (MHD) in its nature. 
We discuss fundamental properties of MHD turbulence. In particular, we discuss the 
generation of compressible MHD waves by Alfvenic turbulence and show that this pro- 
cess is inefficient. This allows us to study the evolution of different types of MHD 
perturbations separately. We describe how to separate MHD fluctuations into 3 dis- 
tinct families - Alfven, slow, and fast modes. We find that the degree of suppression 
of slow and fast modes production by Alfvenic turbulence depends on the strength of 
the mean field. We show that Alfven modes in compressible regime exhibit scalings 
and anisotropy similar to those in incompressible regime. Slow modes passively mimic 
Alfven modes. However, fast modes exhibit isotropy and a scaling similar to that of 
acoustic turbulence both in high and low P plasmas. We show that our findings entail 
important consequences for theories of star formation, cosmic ray propagation, dynam- 
ics of dust, and gamma ray bursts. We anticipate many more applications of the new 
insight to MHD turbulence and expect more revisions of the existing paradigms of 
astrophysical processes as the field matures. 



1 Introduction 

Astrophysics has been providing the major incentive for MHD studies. High 
conductivity of astrophysical fluids makes magnetic fields "frozen in" , and they 
afl'ect fluid motions. The coupled motion of magnetic field and conducting fluid 
is what a researcher has to deal with while studying various astrophysical phe- 
nomena from star formation to gamma ray bursts. 

Turbulence is ubiquitous in astrophysical fluids and it holds the key to many 
astrophysical processes (stability of molecular clouds, heating of the interstellar 
medium, properties of accretion disks, cosmic ray transport etc). Why would we 
expect astrophysical fluids to be turbulent? A fluid of viscosity u gets turbulent 
when the rate of viscous dissipation, which is ^ z^/L^ at the energy injection scale 
L, is much smaller than the energy transfer rate ~ Vl/L, where Vl is the velocity 
dispersion at the scale L. The ratio of the two rates is the Reynolds number 
Re = VlL/v. In general, when Re is larger than 10 — 100 the system becomes 
turbulent. Chaotic structures develop gradually as Re increases, and those with 
Re ~ 10'^ are appreciably less chaotic than those with Re ~ 10*. Observed 
features such as star forming clouds are very chaotic with Re > 10*, which 
ensures that the fluids are turbulent. The measured statistics of fluctuations 
ISM [1,94,16] and Solar wind fluctuations [57] show signatures of the Kolmogorov 
statistics obtained for incompressible unmagnetized turbulent fluid. 
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Kolmogorov theory [40] provides a scaling law for incompressible no7>magnetized 

hydrodynaniic turbulence. This law is true in the statistical sense and it provides 
a relation between the relative velocity vi of fluid elements and their separation 
I, namely, vi <^ l^/^. An equivalent description is to express spectrum E{k) 
as functions of wave number k 1/0- The two descriptions are related by 
kE{k) ~ v"^. The famous Kolmogorov spectrum is E{k) ~ k~^^^. The apphca- 
tions of Kolmogorov theory range from engineering research to meteorology (see 
[75]) but its astrophysical applications are poorly justified. 

Unlike laboratory turbulence astrophysical turbulence is magnetized and 
highly compressible. Then, why do astrophysical fluids show signatures of Kol- 
mogorov statistics? Let us consider incompressible MHD turbulence first. There 
have long been understanding that the MHD turbulence is anisotropic (e.g. [92]). 
A substantial progress has been achieved recently by Goldreich & Sridhar [24] 
(hereafter GS95) who made a prediction regarding relative motions parallel and 
perpendicular to magnetic field B for incompressible MHD turbulence. The GS95 
model envisages a Kolmogorov spectrum of velocity and the scale-dependent 
anisotropy (see below) . These relations have been confirmed numerically (Cho & 
Vishniac [13]; Maron & Goldreich [67]; Cho, Lazarian & Vishniac [10], hereafter 
CLV02b; see also review by Cho, Lazarian, & Vishniac [9], hereafter CLV02a); 
they are in good agreement with observed and inferred astrophysical spectra (see 
CLV02a). A remarkable fact revealed in CLV02b is that fluid motions perpen- 
dicular to B are identical to hydrodynamic motions. This provides an essential 
physical insight into why in some respects MHD turbulence and hydrodynamic 
turbulence are similar, while in other respects they are different. 

However, in most cases compressibility of turbulence is important. For in- 
stance, interstellar medium is highly compressible and star formation requires 
considering supersonic compressible motions (see reviews [99,64,106]). It can be 
shown that assuming that only incompressible turbulence exists in ISM results in 
grossly erroneous conclusions for cosmic ray transport (see review [46]). It may 
be an important question whether the physical pictures in incompressible and 
compressible turbulence are similar. For instance, Cho, Lazarian & Vishniac [11] 
(henceforth CLV02c) reported a new regime of turbulence that takes place in a 
partially ionized gas. In this regime turbulent energy protrudes to small scales 
through a magnetic cascade, while the turbulent velocities are suppressed. How 
will compressibility affect this regime? 

Compressible turbulence is an unsolved problem even in the absence of mag- 
netic flelds. How feasible is it to strive for obtaining universal scaling relations 
for compressible media in view of the fact that no such universality exists for 
compressible hydro turbulence? The difficulty one encounters while studying 
compressible MHD is that MHD turbulence is in general more complicated than 
its hydrodynamic counterpart. In compressible regime, 3 different types of mo- 
tions (Alfven, slow, and fast modes) exist. Alfven modes are incompressible and 
sometimes called shear Alfven modes. The other two modes are compressible 
modes (see §3.2). How do those modes interact? Is it reasonable to talk about 
separate modes in highly non-linear MHD turbulence? These and similar ques- 
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tions dealing with fundamental properties of MHD turbulence we will attempt 
to answer below. 

Thus we must consider a more realistic case of compressible MHD turbulence. 
Literature on the properties of compressible MHD is very rich (see CLV02a). 
Back in 80s Higdon [29] theoretically studied density fluctuations in the inter- 
stellar MHD turbulence. Matthaeus & Brown [68] studied nearly incompress- 
ible MHD at low Mach number and Zank & Matthaeus [105] extended it. In 
an important paper Matthaeus et al. [69] numerically explored anisotropy of 
compressible MHD turbulence. However, those papers do not provide universal 
scalings of the GS95 type. Some hints about effects of compressibility can be 
inferred from Goldreich & Sridhar (GS95) seminal paper. A more focused dis- 
cussion was presented in the Lithwick & Goldreich [61] paper which deals with 
electron density fluctuations in the pressure dominated plasma, i.e. in high (3 
regime (/? = Pgas/Pmag ^ !)• Incompressible regime formally corresponds to 
(3 ^ oo and therefore it is natural to expect that for /3 1 the GS95 picture 
would persist. Lithwick & Goldreich [61] also speculated that for low (3 plasmas 
the GS95 scaling of slow modes may be applicable. A direct study of MHD modes 
in compressible low /3 plasmas is given in Cho & Lazarian [6] (hereafter CL02), 
and more general results applicable for a wide range of (3 and Mach numbers are 
presented in Cho & Lazarian [8] (hereafter CL03). 

The generation of slow and fast modes (i.e. MHD version of sound waves") 
has important astrophysical implications. First, in the presence of damping, 
density and non-Alfvenic magnetic fluctuations are generated only through com- 
pressible daughter waves (i.e. slow and fast waves) generated by Alfven turbu- 
lence. These fluctuations are important for interstellar physics and cosmic ray 
physics. Second, if Alfvenic modes produce a copious amount of compressible 
modes, the whole picture of independent Alfvenic turbulence fails. Therefore, 
ineflacient generation of compressible modes from Alfven turbulence is a neces- 
sary condition for independent Alfvenic cascade. 

In what follows we review observational data on statistics of turbulence, in- 
cluding the velocity data available through spectral line studies §2. In §3 we 
describe our technique for decomposing MHD turbulence into Alfvc^n, slow and 
fast modes. Mode coupling is discussed in §4, while simple theoretical arguments 
about mode scalings are provided in §5. We describe scalings of velocity fluc;- 
tuations in §6 and magnetic and density fluctuations in §7. The new regime of 
turbulence that emerges below the viscous cut-off is briefly discussed in §8. §9 
deals with the applicability of our results and with their significance for the 
theories of star formation, cosmic ray propagation, gamma ray bursts etc. The 
summary is given in §10. 

2 Observational Motivation 

Observations as well as space missions provide data on the statistics of astro- 
physical turbulence. This data suggests that in a wide variety of circumstances 
astrophysical turbulence exhibits power-law spectra consistent with Kolmogorov 
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picture. It would be very naive to think that, in the presence of dynamically im- 
portant magnetic fields, the turbulence may really be Kolmogorov, but it is 
suggestive that, for a wide variety of circumstances, the turbulence should allow 
pretty simple statistical description. This strongly motivates a quest for simple 
relations to describe the apparently complex phenomenon. 

Direct studies of turbulence^ have been done mostly for interstellar medium 
and for the Solar wind. While for the Solar wind in-situ measurements are possi- 
ble, studies of interstellar turbulence require inverse techniques to interpret the 
observational data. 

Attempts to study interstellar turbulence with statistical tools date as far 
back as the 1950s [30,35,78,100] and various directions of research achieved var- 
ious degree of success (see reviews by [36,17,1,43,44,51]). 

2.1 Solar wind 

Solar wind (see review [26]) is a magnetized flow of particles (mostly electrons 
and protons) from the Sim. Studies of the solar wind allow point-wise statistics to 
be measured directly using spacecrafts. These studies are the closest counterpart 
of laboratory measurements. 

The solar wind flows nearly radially away from the Sun, at up to ^700 
km/s. This is much faster than both spacecraft motions and the Alfvcn speed. 
Therefore, the turbulence is "frozen" and the fluctuations at frequency / are 
directly related to fluctuations at the scale k in the direction of the wind, as 
k — 27r//f , where v is the solar wind velocity [32]. 

The solar wind shows f^^^'^ scaling on small scales. The turbulence is strongly 
anisotropic (see [39]) with the ratio of power in motions perpendicular to the 
magnetic field to those parallel to the magnetic field being around 30. The in- 
termittency of the solar wind turbulence is very similar to the intermittency 
observed in hydro dynamic flows [31]. 

2.2 Electron density statistics in the ISM 

Studies of turbulence statistics of ionized media in the interstellar space (see 
[93]) have provided information on the statistics of plasma density at scales 10*- 
10^^ cm. This was based on a clear understanding of processes of scintillations 
and scattering achieved by theorists^ (see [82,27]). A peculiar feature of the 
measured spectrum (see [1]) is the absence of the slope change at the scale at 
which the viscosity by neutrals becomes important. 

Scintillation measurements are the most reliable data in the "big power law" 
plot in Armstrong et al. [1] . However there are intrinsic limitations to the scintil- 
lations technique due to the limited number of sampling directions, its relevance 

^ Indirect studies include the line-velocity relationships [41] where the integrated ve- 
locity profiles are interpreted as the consequence of turbulence. Such studies do not 
provide the statistics of turbulence and their interpretation is very model dependent. 

^ In fact, the theory of scintillations was developed first for the atmospheric applica- 
tions. 
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only to ionized gas at extremely small scales, and the impossibility of getting 
velocity (the most important!) statistics directly. Therefore with the data one 
faces the problem of distinguishing actual turbulence from static density struc- 
tures. Moreover, the scintillation data do not provide the index of turbulence 
directly, but only show that the data arc consistent with Kolmogorov turbu- 
lence. Whether the (3D) index can be -4 instead of -11/3 is still a subject of 
intense debate [29,82]. In physical terms the former corresponds to the superpo- 
sition of random shocks rather than eddies. 

2.3 Velocity and density statistics from spectral lines 

Atoms and molecules in the interstellar space emit radiation at specific wave- 
lengths. A spectral line from atomic hydrogen with Ao=21 cm is particularly 
important in astronomy. Astronomers observe intensity of radiation at different 
wavelengths near Aq for different points on the sky, which results in a 3D data 
cube that consists of two spatial (or, angular) coordinates and one wavelength 
coordinate (i.e. T=T{0i, 02, X — Aq), where T is so-called antenna tempera- 
ture, which measures radiation energy). Using the Doppler shift formula, we can 
convert the wavelength dimension to the velocity dimension. Such spectral line 
data cubes arc unique sources of information on interstellar turbulence. Doppler 
shifts due to supersonic motions contain information on the turbulent velocity 
field which is otherwise difficult to obtain. Moreover, the statistical samples are 
extremely rich and not limited to discrete directions. In addition, line emission 
allows us to study turbulence at large scales, comparable to the scales of star 
formation and energy injection. 

However, the problem of separating velocity and density fluctuations within 
HI data cubes is far from trivial [42,44,49,51]. The analytical description of the 
emissivity statistics of channel maps (velocity slices) in Lazarian & Pogosyan 
[49] (see also [44,51] for reviews) shows that the relative contribution of the 
density and velocity fluctuations depends on the thickness of the velocity slice. 
In particular, the power-law asymptote of the emissivity fluctuations changes 
when the dispersion of the velocity at the scale imder study becomes of the order 
of the velocity slice thickness (the integrated width of the channel map). These 
results are the foundation of the Velocity-Channel Analysis (VGA) technique 
which provides velocity and density statistics using spectral line data cubes. The 
VGA has been successfully tested using data cubes obtained via compressible 
magnetohydrodynamic simulations and has been applied to Galactic and Small 
Magellanic Gloud atomic hydrogen (HI) data [52,49,94,16]. Furthermore, the 
inclusion of absorption effects [50] has increased the power of this technique. 
Finally, the VGA can be applied to different species (GO, Hq, etc.) which should 
further increase its utility in the future. 

Within the present discussion a number of results obtained with the VGA 
are important. First of all, the Small Magellanic Gloud (SMG) HI data exhibit a 
Kolmogorov-type spectrum for velocity and HI density from the smallest resolv- 
able scale of 40 pc to the scale of the SMG itself, i.e. 4 kpc. Similar conclusions 
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can be inferred from the Galactic data [28] for scales of dozens of parsecs, al- 
though the analysis has not been done systematically. Deshpande et al. [16] 
studied absorption of HI on small scales toward Gas A and Gygnus A. Within 
the VGA their results can be interpreted as implying that on scales less than 
1 pc the HI velocity is suppressed by ambipolar drag and the spectrum of density 
fluctuations is shallow P{k) ~ k~'^-^ . Such a spectrum [15] can account for the 
small scale structure of HI observed in absorption. 

2.4 Magnetic field statistics 

Magnetic field statistics arc the most poorly constrained aspect of ISM turbu- 
lence. The polarization of starlight and of the Far-Infrared Radiation (FIR) from 
aligned dust grains is affected by the ambient magnetic fields. Assuming that 
dust grains arc always aligned with their longer axes perpendicular to magnetic 
field (see review [45]), one gets the 2D distribution of the magnetic field direc- 
tions in the sky. Note that the alignment is a highly non-linear process in terms 
of the magnetic field and therefore the magnetic field strength is not available"^. 

The statistics of starlight polarization (see [20]) is rather rich for the Galactic 
plane and it allows to establish the spectrum* E{K) ~ K~^-^, where K is& two 
dimensional wave vector describing the fluctuations over sky patch.''' 

For uniformly sampled turbulence it follows from Lazarian & Shutenkov [53] 
that E{K) ~ K"" for K < Kq and K'^ for K > Kq, where K^^ is the critical 
angular size of fluctuations which is proportional to the ratio of the injection 
energy scale to the size of the turbulent system along the line of sight. For 
Kolmogorov turbulence a = —11/3. 

However, the real observations do not uniformly sample turbulence. Many 
more close stars are present compared to the distant ones. Thus the intermediate 
slops are expected. Indeed, Gho & Lazarian [7] showed through direct simulations 
that the slope obtained in [20] is compatible with the underlying Kolmogorov 
turbulence. At the moment FIR polarimetry does not provide maps that are 
really suitable to study turbulence statistics. This should change soon when 
polarimetry becomes possible using the airborne SOFIA observatory. A better 
understanding of grain alignment (see [45]) is required to interpret the molecular 
cloud magnetic data where some of the dust is known not to be aligned (see [47] 
and references therein). 

Another way to get magnetic field statistics is to use synchrotron emission. 
Both polarization and intensity data can be used. The angular correlation of po- 
larization data [2] shows the power- law spectrum K~^-^ and we believe that the 

^ The exception to this may be the alignment of small grains which can be revealed 
by microwave and UV polarimetry [45] . 

Earlier papers dealt with much poorer samples (see [36]) and they did not reveal 
power-law spectra. 

® This spectrum is obtained by [20] in terms of the expansion over the spherical har- 
monic basis Yim- For sufficiently small areas of the sky analyzed the multipole anal- 
ysis results coincide with the Fourier analysis. 
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interpretation of it is similar to that of starlight polarization. Indeed, Faraday 
depolarization limits the depth of the sampled region. The intensity fluctuations 
were studied in [53] with rather poor initial data and the results were inconclu- 
sive. Cho & Lazarian [7] interpreted the fluctuations of synchrotron emissivity 
[21,22] in terms of turbulence with Kolmogorov spectrum. 

3 Numerical Approach 

3.1 Helmholtz decomposition for hydrodynamic turbulence 

To get an insight of the turbulence cascade we have attempted a decomposition 
of the MHD turbulent flow into Alfven, slow and fast modes (see CL02, CL03). 
Our numerical method is similar to the technique utilizing the ^^Helmholtf de- 
composition in hydrodynamics. 

Our method is different from Lighthill's theory [60] of Jar field acoustic 
wave generation from homogeneous turbulence. Literature on the application 
of Lighthill's approximation to astrophysical problems is rich. Astrophysical flu- 
ids are stratified (by gravity) and/or magnetized. Therefore, Lighthill's theory 
requires modifications for astrophysical fluids. Stein [96] extended Lighthill's 
theory to stratified astrophysical fiuids in gravitational field. Subsequent papers 
(e.g. [23]) further discussed about generation of acoustic waves in (Solar) con- 
vection zone. On the other hand, Musielak fc Rosner [79] constructed a model 
for weak magnetic field convection zone and Musielak, Rosner, & Ulmschneider 
[80] discussed about wave generation in an inactive flux tube. Lee [58] explored 
wave generation in sunspots, which are the most strongly magnetized on the 
surface of the Sun. All these approaches are to calculate far-field acoustic flux. 

Moyal [76] introduced a method that decomposes velocity fleld in Fourier 
space, which is equivalent to Helmholtz's decomposition of a vector fleld: V = 

-I- Vc, where Vg is divergence- free (V • Vs=0) field and Vc is curl-free 
(V X Vc=0) field. Note that Vg represents incompressible or solenoidal part 
and Vc compressible or dilatational one. In Fourier space, solenodal and di- 
latational components have simple geometrical meanings: Vg is the component 
perpendicular to the wave vector k and Vg parallel to k. 

The first (published) numerical simulations of compressible hydrodynamic 
turbulence were performed by Feiereisen, et al. [19]. They studied subsonic (sonic 
Mach numbers, M^, up to 0.32) homogeneous shear fiows with 64'^ grid points. 
Passot & Pouquet [83] carried out two-dimensional isotropic homogeneous com- 
pressible decaying turbulence with 256^ grid points. They showed that properties 
of turbulence at low initial Mach numbers [Mg < 0.3) is significantly different 
from those of higher Mach number ones. Passot, Pouquet, & Woodward [84] sim- 
ulated two-dimensional isotropic decaying tm-bulenc;e with initial Mach numbers 
up to 4. They provided conjecture for the three-dimensional case and discussed 
implications of their work on astrophysical fluids in the interstellar medium. 
Subsequent simulations [37,38,95,89,59] addressed various issues of compressible 
turbulence. Recent high resolution three-dimensional simulations include Porter, 
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Pouquet, & Woodward [86], Porter, Woodward, & Pouquet [88], and Porter, 
Pouquct, & Woodward [87]. 

The energy spectra of compressible hydrodynamic turbulence are still uncer- 
tain. For spectrum of solenoidal components, a Kolmogorov-type dimensional 
analysis leads to 

Esolerundalik) CC k-^'l^ikLT^" (1) 

([34]; see also [84]). However, Moiseev et al. [74] obtained slightly different results 

Esolenoidal{k)^k-''/\kL)^^^K (2) 

When, Ms 0, both results give Kolmogorov spectrum. 

The energy spectrum of compressible components is more uncertain. For 
example, Zakharov & Sagdeev [104] derived scalings for compressible modes: 

Erad{k) (X (3) 

where the subscript rad denotes compressible components (i.e. radial compo- 
nents in Fourier space). On the other hand, Bataille & Zhou [3] and Bertoglio, 
Bataille, & Marion [4] obtained that the spectral index (slope) is a function of 
Mach number, Mg. When, Mach number is of order unity, their results give a 
Kolmogorov spectrum. Recent numerical simulations [88] with up to 1024^ grid 
points show Kolmogorov's k~^^^ spectra both for Esoienoidai{k) and Erad{k)- 

The generation of compressible components from incompressible initial tur- 
bulence is also an unresolved issue. Closure calculation by Bataille & Zhou [3] 
and Bertoglio et al. [4] predicts that x = (V;L)/(Koieno.dai) ~ ^1 Numeri- 
cal calculations of decaying turbulence with initial Mach number of order unity 
[86,87] show that x ~ 0.1. 

3.2 MHD mode decomposition 

Three types of waves exist (Alfven, slow and fast) in compressible magnetized 
plasma. In this section, we describe how to separate different MHD modes. 

In the presence of magnetic field B, the momentum equation has an addi- 
tional term, (V x B) x B (divided by l/47r). This is the so-called J x B term, 
which can be re-written as the sum of the magnetic tension term, B • VB, and 
magnetic pressure term, V^B: 

(V X B) X B = B • VB - V2B/2. (4) 

In addition, when magnetic Reynolds number, Rm = VL/t], where rj is magnetic 
diffusivity, is large, magnetic field lines move together with fluid elements, which 
is sometime called that magnetic fields are frozen-in. 

In some sense, the magnetic field lines are like elastic bands moving together 
with fluid elements in that they have tension. However, they are different from 
rubber bands in that they are repulsive each other, which is the nature of mag- 
netic pressure. 
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Because of tension and pressure, the nature of MHD waves is much more 

complicated than their hydrodynamic counterpart - sound wave. This is because 
we need to consider 3 different restoring forces - magnetic tension, magnetic 
pressure, and gas pressure. For Alfven waves, magnetic tension is the only the 
restoring force (Fig. 1(a)). For slow and fast waves, all 3 restoring forces are 
important. For Slow waves, magnetic and gas pressure are out of phase and, for 
fast modes, they are in phase (Fig. 2). 

The slow, fast, and Alfven bases that denote the direction of displacement 
vectors for each mode arc given by 



oc (-1 + a- Vi:')fc||k|j 
if cx (-1 + a + y/D)k\\k\\ 
|a oc cx k_L X k|| , 



(1- 
(1 



a - 
- a - 



(5) 
(6) 
(7) 



where D — {1 + a)^ — 4acos(?, a = a^/Vj = P{j/2), 9 is the angle between k 
and Bq, and ip is the azimuthal basis in the spherical polar coordinate system 
(see Appendix). These are equivalent to the expression in CL02: 



oc fc||k|| + 



if « 



l-\/D-/3/2 
l + VD + f3/2 



k± 



A;_Lkj 



1 - J* + P/2 
1 + VD-(3/2 



k\i 



fc||k|| + fc_Lk_L. 



(8) 
(9) 



(Note that 7 = 1 for isothermal case.) 

We can obtain slow and fast velocity by projecting velocity Fourier compo- 
nent Vk into is and respectively. In Appendix, we also discuss how to separate 
slow and fast magnetic modes. We obtain energy spectra using this projection 
method. 



3.3 Numerical method 

We use a hybrid essentially non-oscillatory (ENO) scheme to solve the ideal 
isothermal MHD equations. When variables arc sufficiently smooth, we use the 
3rd-order Weighted ENO scheme [33] without characteristic mode decomposi- 
tion. When the opposite is true, we use the 3rd-order Convex ENO scheme [62] . 
Combined with a three-stage Runge-Kutta method for time integration, our 
scheme gives third order accuracy in space and time. We solve the ideal MHD 
equations in a periodic box: 

dp/dt + V • (pv) = 0, (10) 
dv/dt -I- V • Vv -I- p- V(aV) - (V x B) x B/Airp = f , (11) 
dB/dt - V X (v X B) = 0, (12) 

with V-B = and an isothermal equation of state. Here f is a random large-scale 
driving force, p is density, v is the velocity, and B is magnetic field. The rms 



10 Cho & Lazarian 




Real space Fourier space 

(c) 

Fig. 1. (a) Alfven wave and its direction of oscillation/displacement. Tlie direction of 
displacement is perpendicular to both Bq and k. Only magnetic tension is the restor- 
ing force, (b) Slow wave and its direction of osciUation/displacement. The direction 
of displacement is between —6 and Bo- (c) Slow wave and its direction of oscilla- 
tion/displacement. The direction of displacement is between k and kj_. 
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SLOW WAVE ^ % ^ FAST WAVE 




Fig. 2. Slow and fast waves in real space (Bo — k plane). We show the directions of 
displacement vectors for a slow wave {left panel) and a fast wave {right panel) . Note that 
lies between d and Bo (= k||) and ^/ between k and kx- Again, 6 is perpendicular 
to k and parallel to the wave front. Note also that, for the fast wave, for example, 
density (inferred by the directions of the displacement vectors) becomes higher where 
field lines are closer, resulting in a strong restoring force, which is why fast waves are 
faster than slow waves. From CLV02a. 

velocity SV is maintained to be unity, so that v can be viewed as the velocity 
measured in units of the r.m.s. velocity of the system and 'B/^/A^Tp as the Alfven 
speed in the same units. The time t is in units of the large eddy turnover time 
(~ L/SV) and the length in units of L, the scale of the energy injection. The 
magnetic field consists of the uniform background field and a fluctuating field: 
B = Bo + b. 

For mode couphng studies (Fig. 4), we do not drive turbulence. For scahng 
studies, we drive turbulence solenoidally in Fourier space and use 216^ points 
and pq = 1. The average rms velocity in statistically stationary state is 6V ~ 0.7. 

For our calculations we assume that Bq / y^inp ~ SB/y/iirp ~ SV. In this 
case, the sound speed is the controlling parameter and basically two regimes 
can exist: supersonic and subsonic. Note that supersonic means low-beta and 
subsonic means high-beta. When supersonic, we consider mildly supersonic (or, 
mildly low-/?) and highly supersonic (or, very \ow-(3). 
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Fig. 3. Separation method. We separate Alfven, slow, and fast modes in Fourier space 
by projecting the velocity Fourier component Vk onto bases ^a, ^s, and respectively. 
Note that = —>f>- Slow basis and fast basis ^/ lie in the plane defined by Bo and 
k. Slow basis lies between —6 and ky. Fast basis ^/ lies between k and kx. Prom 
CL03. 



4 Mode Coupling: Theory and Simulations 

As mentioned above, the coupling of compressible and incompressible modes is 
crucial. If Alfvenic modes produce a copious amount of compressible modes, the 
whole picture of independent Alfvenic turbulence fails. 

The generation of compressible motions (i.e. radial components in Fourier 
space) from Alfvenic turbulence is a measure of mode coupling. How much energy 
in compressible motions is drained from Alfvenic cascade? According to the 
closure calculations ([4]; see also [105]), the energy in compressible modes in 
hydrodynamic turbulence scales as ^ if Af, < 1. We may conjecture that 
this relation can be extended to MHD turbulence if, instead of M^, we use 
~ {SV)\/{a'^ + V|). (Hereinafter, we define Va ^ Bo/^inp.) However, as the 
Alfven modes are anisotropic, this formula may require an additional factor. 
The compressible modes are generated inside so-called Goldreich-Sridhar cone, 
which takes up ^ {SV)a/Va of the wave vector space. The ratio of compressible 
to Alfvenic energy inside this cone is the ratio given above. If the generated 
fast modes become isotropic (see below), the diffusion or, "isotropization" of 
fast wave energy in the wave vector space increase their energy by a factor of 
~ VA/iSV)A- This results in 



2 

rad 



ml 



m\ Va 



(13) 



where (SV)"^^^ and {SV)a are energy of compressible^ and Alfven modes, re- 
spectively. Eq. (13) suggests that the drain of energy from Alfvenic cascade is 
marginal when the amplitudes of perturbations are weak, i.e. {6V)a *C Va- 

® It is possible to show that the compressible modes inside the Goldreich-Sridhar cone 
are basically fast modes. 
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Fig. 4(a) shows that generation of slow and fast modes (the dotted hne) from 
Alfven modes (the sohd hnc) is marginal. The result shown in the figure is for 
Ms = 1.6 at i = 0. We repeated similar simulation for different Mach numbers 
and plasma /3's and measured the ratios of energy in compressible modes to that 
in Alfven modes. The results arc shown in Fig. 4(c) and (d). Fig. 4(c) suggest 
that the generation of compressible motions follows equation (13). Fast modes 
also follow a similar scaling, although the scatter is a bit larger. Fig. 4(e) demon- 
strates that fast modes arc initially generated anisotropically, which supports our 
theoretical consideration above. Fast modes becomes isotropic later (Fig. 4(f)). 
Fig. 4(d) shows that generation of slow modes follows {SV)1/{6V)\ oc {6V)a/Va 
for low (3 cases (pluses in the figure). But, the scaling is not clear for high /3 cases 
(diamonds in the figure). 

Fig. 4(b) shows that dynamics of Alfven modes is not aff'ected by slow modes. 
The solid line in the figure is the energy in Alfven modes when we start the decay 
simulation with Alfven modes only. The dotted line is the Alfven energy when 
we start the simulation with all modes. This result confirms that Alfven modes 
cascade is almost independent of slow and fast modes. In this sense, coupling 
between Alfven and other modes is weak. 

5 Quest for Scaling Relations 

5.1 Scaling of incompressible MHD turbulence 

As mentioned in §1, Goldreich & Sridhar (GS95) made a prediction regarding rel- 
ative motions parallel and perpendicular to magnetic field B for incompressible 
MHD turbulence. Here, we reconstruct GS95 model from different perspectives. 

An important observation that leads to understanding of the GS95 scaling 
is that magnetic field cannot prevent mixing motions of magnetic field lines 
if the motions are perpendicular to magnetic field. Those motions will cause, 
however, waves that will propagate along magnetic field lines. If that is the case, 
the time scale of wave-like motion, i.e. ^ /||/V/i, where Z|| is the characteristic 
size of the perturbation and Va is the local Alfven speed, will be equal to the 
hydrodynamic time-scale, l±/vi. The mixing motions are hydrodynamic-like and 

1/3 

therefore obey Kolmogorov scaling vi oc . Equating the two relations above, 
we obtain a critically balance condition 

ll\/VA^l±/vk (oT k\\VA ^ k±Vk). (14) 

If conservation of energy in the turbulent cascade applies locally in phase space 
then the energy cascade rate (vf /teas) is constant): (vf)/(l±/vi) = constant. 
Combining this with the critical balance condition we obtain 

Z|| oc lY^ (15) 

(or fcji oc k'^^ in terms of wavevectors) and and a Kolmogorov-like spectrum for 
perpendicular motions 

VI oc or, E(k) oc k~^^^, (16) 




Fig. 4. Generation of compressible MHD modes, (a) When Alfvenic turbulence decays, 
generation of slow and fast modes is marginal. 144"'. Initially, /3 (ratio of gas to magnetic 
pressure, Pg/Pmag) = 0.2 and Ms (sonic Mach number) ~ 1.6. (b) Comparison of decay 
rates. Decay of Alfven modes is not much affected by other (slow and fast) modes. We 
use 144^ grid points. Initially, /3 = 0.2 and Ma ~ 1.6 for the solid line and Ma ~ 2.3 
for the dotted line, (c) Square of the r.m.s. velocity of the compressible modes. We use 
144^ grid points. Only Alfven modes are allowed as the initial condition. "Pluses" are 
for low P cases (0.02 < (3 < 0.4). "Diamonds" are for high /? cases (1 < /3 < 20). Fast 
modes follow a similar scaling, (d) Square of the r.m.s. velocity of the slow modes. See 
(c) for explanation, (e) Generation of fast modes. Snapshot is taken at t=0.06 from a 
simulation (with 144'' grid points) that started off with Alfven modes only. Initially, 
(3 = 0.2 and Ms (sonic Mach number) ^ 1.6. (f) Generation of fast modes. Snapshot 
is taken at t~2. See (e) for explanation. 
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which is not surprising since perpendicular motions are hydrodynamic. If we 

interpret Z[| as the eddy size in the direction of the locaF field and l± as that in the 
perpendicular direction, the relation in equation (15) implies that smaller eddies 
are more elongated (see Fig. 5 for illustration of scale-dependent anisotropy). 

Using Matthaeus et al. [71] result, we can re-derive GS95 model. Afatthaeus 
et al. [71] showed that the anisotropy of low frequency MHD turbulence scales 
linearly with the ratio of perturbed and total magnetic field strength b/B {= 
b/ib"^ + Bq)^/"^). This scaling relation has simple geometric meaning: perpendic- 
ular size of a large scale eddy is similar to its parallel size times b/B, which 
is is determined by magnetic field line wandering. Although their analysis was 
based on comparing the strength of a uniform background field and the mag- 
netic perturbations on all scales, we can reinterpret this result by assuming that 
the strength of random magnetic field at a scale I is 6;, and that the back- 
ground field is the sum of all contributions from larger scales. Then Matthaeus 
et al.'s result becomes a prediction that the anisotropy (fcy/fcj,) is proportional 
to {bi/B). We can take the total magnetic field strength B ^ constant as long 
as the background field is stronger than the perturbations on all scales. Since 
bi ~ (kE{k)Y/'^ ~ A; _|_^^'^, we obtain an anisotropy (A;||/fc_L) proportional to k']^^^ , 

and oc k]^^ . In this interpretation, smaller eddies are more elongated because 
they have a smaller bi/B ratio. 

5.2 Compressible scalings: theoretical considerations 

In §4, we showed that Alfven modes are independent of other modes. Therefore, 
we expect GS95 scalings for Alfven modes even for supersonic turbulence. 

When Alfven cascade evolves on its own, it is natural to assume that slow 
modes passively follow the Alfven cascade and exhibit the same scaling. Indeed, 
slow modes in high (3 plasmas are similar to the pseudo- Alfven modes in incom- 
pressible regime (see GS95; [61]). The latter modes do follow the GS95 scaling. 
In low (3 plasmas, motions of slow modes are density perturbations propagating 
with the sound velocity a parallel to the mean magnetic field (see equation (76)). 
In magnetically dominated environments (/3 <C 1), a <C Va and the gaseous per- 
turbations are essentially static. Therefore the magnetic field mixing motions 
are expected to mix density perturbations as if they were passive scalar. It is 
known that the passive scalar shows the same scaling as the velocity field of the 
inducing turbulent motions. Thus the slow waves are expected to demonstrate 
GS95 scafings (see CL02). 

The fast waves in low /3 regime propagate at Va irrespectively of the magnetic 
field direction. In high /3 regime, the properties of fast modes are similar, but 
the propagation speed is the sound speed a. Thus the mixing motions induced 
by Alfven waves should affect the fast wave cascade only marginally. The latter 
cascade is expected to be analogous to the acoustic wave cascade and be isotropic. 

'' The concept of local is crucial. The GS95 scalings are obtained only in the local frame 
of magnetic field, as this is the frame where magnetic field are allowed to be mixed 
without being opposed by magnetic tension. 
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6 Velocity scaling 

6.1 Illustration of eddy structures 

Fig. 5 and Fig. 6 show the shapes of eddies of different sizes. For Alfven mode 
eddies (Fig. 5), left 3 panels show an increased anisotropy as we move from the 
top (large eddies) to the bottom (small eddies). The horizontal axes of the left 
panels are parallel to Bq. Structures in the perpendicular plane (right panels) 
do not show a systematic elongation. However, Fig. 6 shows that velocity of 
fast modes exhibit isotropy. Data are from a simulation with 216^ grid points, 
Ms = 2.3, and (3 = 0.2. 

6.2 Alfven modes in compressible MHD 

If Alfven cascade evolves on its own, it is natural to assume that slow modes 
exhibit the GS95 scaling. Indeed, slow modes in pressure dominated environ- 
ment (high j3 plasmas) are similar to the pseudo- Alfven modes in incompressible 
regime (see GS95; [61]). The latter modes do follow the GS95 scaling. In mag- 
netically dominated environments (low (3 plasmas), slow modes arc density per- 
turbations propagating with the sound velocity a parallel to the mean magnetic 
field (see equation (76)). Those perturbations are essentially static for a <C Va- 
Therefore Alfvenic turbulence is expected to mix density perturbations as if they 
were passive scalar. This also induces GS95 spectrum. 

Fig. 7(a), (c), and (e) show that the spectra of Alfven waves follow a Kol- 
mogorov spectrum: 

Spectrum of Alfven Waves: E^{k) oc kj^^^, (17) 

regardless of plasma (3 or sonic Mach number Mg. 

In Fig. 7(b), (d), and (f), we plot contours of equal second-order structure 
function for velocity (5F2(r) =< |v(x -|- r) — v(x)p >avg. over x) obtained in 
local coordinate systems in which the parallel axis is aligned with the local 
mean field (see [13];[67]; CLV02b). The SF2 along the axis perpendicular to the 
local mean magnetic field follows a scaling compatible with r^^^. The SF2 along 
the axis parallel to the local mean field follows steeper scaling. The results 
are compatible with the GS95 model, 

Anisotropy of Alfven Waves: ru oc r'^^, or fc|| oc fc^/^, (18) 

where r|| and r± are the semi-major axis and semi-minor axis of eddies, respec- 
tively [13]. When we interpret that contours represent eddy shapes, the above 
scaling means that smaller contours are more elongated. 

6.3 Slow modes in compressible MHD 

The incompressible limit of slow waves is pseudo-Alfven waves. Goldreich & 
Sridhar [25] argued that the pseudo-Alfven waves are slaved to the shear- Alfven 
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Fig. 5. Anisotropy as a function of scale. Alfven mode velocity show scale-dependent 
anisotropy. Lighter tones are for larger |v|. Only part of the data cube is shown. Mag- 
netic field show similar behaviors. Large scale eddies are obtained from the Fourier 
components with 1 < fc < 3. Medium scale eddies are obtained from the Fourier com- 
ponents with 3 < fc < 9. Small scale eddies are obtained from the Fourier components 
with 9 < fc < 27. 



(i.e. ordinary Alfven) waves, which means that pseudo- Alfven modes do not cas- 
cade energy for themselves. Lithwick & Goldreich [61] made similar theoretical 
arguments for high (3 plasmas and conjectured similar behaviors of slow modes 
in low (3 plasmas. We confirmed that similar arguments are also applicable to 
slow waves in low (3 plasmas (CL02). Indeed, energy spectra in Fig. 8(a) and (c) 
are consistent with: 

Spectrum of Slow Modes: E^{k) cx k^^^^^. (19) 
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However, the kinetic energy spectrum for sfow modes in Fig. 8(e) does not show 
the Kolmogorov slope. The sfope is close to —2, which is suggestive of shock 
formation. At this moment, it is not clear whether or not the —2 slope is the 
true slope. In other words, the observed —2 slope might be due to the limited 
numerical resolution. Runs with higher numerical resolution should give the 
definite answer. 

In Fig. 8(b), (d), and (f), contours of equal second-order velocity struc- 
ture function {SF2), representing eddy shapes, show scale-dependent anisotropy: 
smaller eddies are more elongated. The results are compatible with the GS95 
model 

Anisotropy of Slow Modes: fc|| oc fc^^'^, or r|| cx r^/"^, (20) 
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Fig. 7. Alfvon modes. (a)&(b): Ms - 0.35 (/3 = 4). (c)&(d): ~ 2.3 {(3 = 0.2). 
(e)&(f): Ms ~ 7 (/3 = 0.02). Spectra are compatible with Kolmogorov. Contours 
(or, eddy shapes) show scale-dependent anisotropy: smaller eddies are more elongated. 
Overall, the scalings are consistent with Goldreich & Sridhar scalings. 



where ry and r± are the semi-major axis and semi-minor axis of eddies, respec- 
tively. 




Fig. 8. Slow modes. (a)&(b): Ms ^ 0.35 (/3 = 4). (c)&(d): Ms ~ 2.3 (/3 = 0.2). 
(e)&(f): Ms ~ 7 (/3 = 0.02). In (a) and (c), spectra are compatible witli Kolmogorov. In 
(e), spectra is uncertain. Contours (or, eddy shapes) show scale-dependent anisotropy: 
smaller eddies are more elongated. Overall, the scalings are consistent with Goldreich 
& Sridhar scalings (except spectra in (e)). 



6.4 Fast modes in compressible MHD 

Fig. 9(b), (d), and (f) show fast modes are isotropic. The resonance conditions 
for the interacting fast waves are Wi + u)2 = and ki + k2 = ka. Since 
a; oc fc for the fast modes, the resonance conditions can be met only when all 
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three k vectors are collinear. This means that the direction of energy cascade is 

radial in Fourier space. This is very similar to acoustic turbulence, turbulence 
caused by interacting sound waves [103,104,63]. Zakharov & Sagdeev [104] found 
E{k) oc k~^/'^. However, there is debate about the exact scaUng of acoustic 
turbulence. Here we cautiously claim that our numerical results are compatible 
with the Zakharov & Sagdeev scaling: 

SpecArum of Fast Modes: E^' (k) k~^^'^ . (21) 

The eddies are isotropic (see also Fig. 6). 



7 Magnetic Field and Density Scalings 

We expect that isotropy/anisotropy of magnetic field is similar to that of velocity 
(see CL03). However, anisotropy of density shows different behavior. Density 
shows anisotropy for the high /? case. But, for low /? cases, density shows more 
or less isotropic structures. We suspect that shock formation is responsible for 
isotropization of density. 

To estimate the r.m.s. fluctuations, we use the following linearized continuity 
and induction equations: 

\pk\ = {poVk/c)\k-i\, (22) 

\bk\ = {BoVk/c)\Bo X i\, (23) 

where c denotes velocity of slow or fast waves (equation (67)). Prom this, we 
obtain the r.m.s. fluctuations 

{Sp/po)s =iSVU\^-i/c,\), (24) 

iSp/po)f ={5V)f{\k-if/cf\), (25) 

idB/Bo)s = {5V)s{\Bo X i/cs\), (26) 

{SB/Bo)f^{SV)f{\Boxij/cf\), (27) 

where angled brackets denote a proper Fourier space average. Generation of slow 
and fast modes velocity {{5V)s and {5V)f) depends on driving force. Therefore, 
we may simply assume that 

{SV)A^{SV)sr^{6V)f, (28) 

where wc ignore constants of ordc;r unity. However, when we consider mostly 
incompressible driving, the generation fast modes may follow equation (13). 
In this case, the amplitude of fast mode velocity is reduced by a factor of 

1/2 



Vi+a/(6VlA 
WIT Va 



{6V)a ~ {5V)s 



{5V)f. (29) 



When we assume {5V)a ~ Bq, equation (29) reduces to equation (28) in low (3 
plasmas. 




Fig. 9. Fast modes. (a)&(b): M„ ~ 0.35 (/3 = 4). (c)&:(d): ~ 2.3 {13 = 0.2). (e)&(f): 
Ms ~ 7 (/3 = 0.02). Spectra are compatible witli either k~^^^ or Kolmogorov. Contours 
(or, eddy shapes) show isotropy. Overall, the scalings are consistent with those of 
acoustic turbulence. 



7.1 L0W-/3 case 

In this limit, Cg ^ acosO and c/ ~ Va- Using equations (76) and (77), we obtain 



{5p/po)s ~ {SV)s{\ cos0/cs\) ~ {SV)s/a, 
{5p/po)f = {6V)f{\sme/cf\) ~ {6V)f/VA, 



(30) 
(31) 
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{SB/Bo)s = {SVU\acosesme/cs\) ~ a{5V)s/a, (32) 

{SB/Bo)f = {SV)f{\l/cf\) ~ {SV)f/VA, (33) 

where we ignore cos ^'s or sin ^'s. 

When we assume {SV)a ^ {SV)s ^ {SV)f ^ Bq, we get 

{dp/po)s (34) 

i6p/po)f =v^M„ (35) 

{SB /Bo), = (36) 

{6B/Bo)f = ^/pMs (37) 



Therefore, in low /? plasmas, slow modes give rise to most of density fluctuations 
(CL02). On the other hand, magnetic fluctuation by slow modes is smaller than 
that by fast modes by a factor of ^/^. 



7.2 High-/3 case 

In this limit, Cg ~ Va cos^ and c/ ~ a. Using equations (78) and (79), we obtain 

(Sp/po). ~ {6V),{\cos9sin0/iac,)\) 

~ {VA/a){dV)s/a, (38) 

i5p/po)f ={SV)f{\l/cf\)^{SV)f/a, (39) 

iSB/Bo)s = iSV)s{\ cos9/cs\) ~ iSV),/VA, (40) 

{SB/Bo)f = {SV) f{\ sine /cf\) ~ {SV)f/a, (41) 

where we ignore cost's or sin6''s. 

Let us just assume that {SV)a ~ (^^)s ~ .Bo ~ -^s~^(^^)/ (cf- equation 
(29)). Then we have 

{Sp/po), ^M./^^Ml (42) 

(5p/po)/ ~M2, (43) 

{5B/Bo)s = 0(1), (44) 

{5B/Bo)f = Ml (45) 

The density fluctuation associated with slow modes is ~ M^, when {SV)s ~ 
{SV)a ^ Va- This is consistent with Zank & Matthaeus [105]. The ratio of {5p)s 
to {5p) f is of order unity. Therefore, both slow and fast modes give rise to similar 
amount of density fluctuations. Note that this argument is of order-of-magnitude 
in nature. In fact, in our simulations for the high /? case, the r.m.s. density 
fluctuation by slow modes is about twice as large as that by fast modes. When 
we use equation (28), we have a different result: {Sp)s ^ {VA/a){5p)f < {Sp)f. 
It is obvious that slow modes dominate magnetic fluctuations: {6B)s > {dB)f 
for both equations (28) and (29). 
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8 Slowly Evolving Fluctuations Below Viscous Cutoff 

In hydrodynamic turbulence viscosity sets a minimal scale for motion, with an 
exponential suppression of motion on smaller scales. Below the viscous cutoff 
the kinetic energy contained in a wavcnumbcr band is dissipated at that scale, 
instead of being transferred to smaller scales. This means the end of the hydro- 
dynamic cascade, but in MHD turbulence this is not the end of magnetic struc- 
ture evolution. For viscosity much larger than resistivity, ^ 3> ?7, there will be a 
broad range of scales where viscosity is important but resistivity is not. On these 
scales magnetic field structures will be created by the shear from non-damped 
turbulent motions, which amounts essentially to the shear from the smallest 
undamped scales. The created magnetic structures would evolve through gener- 
ating small scale motions. As a result, we expect a power-law tail in the energy 
distribution, rather than an exponential cutoff. This completely new regime for 
MHD turbulence was reported in CLV02c. Further research showed that there 
is a smooth connection between this regime and small scale turbulent dynamo 
in high Prandtl number fluids (see [90]). 

CLV02c explored this regime numerically with a grid of 384'^ and a physical 
viscosity for velocity damping. The kinetic Reynolds number was around 100. 
We achieved a very small magnetic diffusivity by the use of hyper-diffusion. The 
result is presented in Fig. 10a. A theoretical model for this new regime and its 
consequences for stochastic reconnection [54] can be found in Lazarian, Vishniac, 
& Cho [55]. It explains the spectrum E{k) ^ as a cascade of magnetic energy 
to small scales under the influence of shear at the marginally damped scales. The 
mechanism is based on the solenoidal motions and therefore the compressibility 
should not alter the physics of this regime of turbulence. 




Fig. 10. Viscously damped regime (viscosity > magnetic diffusivity). Due to large 
viscosity, velocity damps after k ^ 10. (a) Left: Incompressible case with 384^ grid 
points. Magnetic spectra show a shallower slope {Eb{k) oc k~^) below the velocity 
damping scale. From CLV02c. (b) Middle: Compressible case with 216^ grid points 
(high /3). Magnetic and density spectra show structures below the velocity damping 
scale at ~ 10. The spectral slope is poorly defined because it was relatively hard to 
achieve very small magnetic diffusivity in the compressible run. From CL03. (c) Right: 
Compressible case with low /?. Density structures are vividly enhanced. 

CL03 showed that the new regime of turbulence is also valid for compress- 
ible MHD. We use the same physical viscosity as in incompressible case (see 
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CLV02c). We rely on numerical diffusion, which is much smaller than physical 

viscosity, for magnetic field. The incrtial range is much smaller due to numerical 
reasons, but it is clear that the new regime of MHD turbulence persists. The 
magnetic fluctuations, however, compress the gas and thus cause fluctuations in 
density. The amplitude of density perturbations is higher in low (3 plasma (see 
Fig. 10c). This is a new (although expected) phenomenon compared to our ear- 
lier incompressible calculations. These density fluctuations may have important 
consequences for the small scale structure of the ISM. 

9 Discussion 

9.1 Range of applicability 

One may argue that for the first time ever we have universal scaling relations 
that describe turbulence over a wide range of plasma /? and Mach numbers. This 
remarkable fact entails a lot of astrophysical implications. 

Is this consistent with the observational data discussed above? 
Yes, we would claim that the observed spectra that are consistent with the 
Kolmogorov scaling arise naturally. Indeed, Alfven and slow modes exhibit GS95 
type scaling. Most of the energy in those modes correspond to perturbations 
perpendicular to magnetic field. The scaling for the perpendicular motions is 
Kolmogorov- like, i.e. E{k) ~ kj_^^^. However, this fact does not imply that it 
is OK to use Kolmogorov scalings to solve astrophysical problems (see §9.2). 
It can be a serious mistake to disregard scale-dependent anisotropy of MHD 
turbulence. Fast modes are isotropic and have a bit different scaling. A more 
careful analysis of the observational data is necessary to detect the signature of 
fast modes. 

Is super- Alfvenic turbulence different? 
In the paper above we considered the cases in which the Alfven speed associ- 
ated with the mean magnetic field is slightly faster than the r.m.s. fluid velocity. 
This regime is called "sub- Alfvenic" regime. If initially the turbulent energy is 
larger than magnetic energy, we are in the regime of so-called "super- Alfvenic" 
turbulence. In this regime the growth of the magnetic field is expected through 
so called "turbulent dynamo" (see [12]; [77]; CLV02a). The magnetic energy at 
scale / in this regime grows exponentially with the characteristic rate of the 
eddy turnover time. Thus we expect to reach equipartition between magnetic 
and kinetic energies at the energy injection scale. At smaller scales the turbu- 
lence becomes sub- Alfvenic and our earlier considerations should be applicable. 
Although the decomposition of MHD turbulence described above does not work 
for irregular magnetic field characterizing super-Alfvenic turbulence, CL03 re- 
sults are suggestive that our considerations about Alfven, slow and fast modes 
are applicable to this regime. 
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9.2 Compressible MHD turbulence and star formation 

How fast does MHD turbulence decay? 

This question has fundamental imphcations for star form,ation (sec [66]). In- 
deed, it was thought originally that magnetic fields would prevent turbulence 
from fast decay. Later (see [65,97]; and review [99]) this was reported not to be 
true. However, fast decay was erroneously associated with the coupling between 
compressible and incompressible modes. The idea was that incompressible mo- 
tions quickly transfer their energy to the compressible modes, which get damped 
fast by direct dissipation (presumably through shock formation). 

Our calculations support the relation in eq. (13). According to it the coupling 
of Alfven and compressible motions is important only at the energy injection 
scales where SVi ^ Va- As the turbulence evolves the perturbations become 
smaller and the coupling less efficient. Typically for numerical simulations the 
inertial range is rather small and this could explain why marginal coupling of 
modes was not noticed. 

Our results show that MHD turbulence damping does not depend on whether 
the fluid is compressible or not. The incompressible motions damp also within 
one eddy turnover time. This is the consequence of the fact that within the strong 
turbulence^ mixing motions perpendicular to magnetic field are hydrodynamic 
to high order (CLV02b) and the cascade of energy induced by those motions is 
similar to the hydrodynamic one, i.e. energy cascade happens within an eddy 
turnover time. 

Is the decay always fast for compressible MHD turbulence ? 
This issue does need further investigation. However, our preliminary answer 
to this question is "no". Indeed, incompressible MHD computations (see [67]; 
CLV02b) show that the rate of turbulence decay depends on the degree of tur- 
bulence imbalance^, i.e. the difference in the energy fluxes moving in opposite 
directions. The strongly imbalanced incompressible turbulence was shown to 
persist longer than its balanced counterpart. This enabled CLV02b to speculate 
that this may enable energy transfer between clouds and may explain the ob- 
served turbulent linewidths of GMCs without evident star formation. Our results 
above show a marginal coupling of compressible and incompressible modes. This 
is suggestive that the results obtained in incompressible simulations are appli- 
cable to compressible environments if amplitudes of perturbations are not large. 
The complication arises from the existence of the parametric instability [14] that 
happens as the density perturbations reflect Alfven waves and grow in amplitude. 
This instability eventually controls the degree of imbalance that is achievable. 
However, the growth rate of the instability is substantially slower than the Alfven 
wave oscillation rate. Therefore, if we take into account that interstellar sources 

* For a formal definition of strong, weak and intermediate turbulence see Goldreich & 
Sridhar [25] and CLV02a, but here we just mention in passing that in most astro- 
physically important cases the MHD turbulence is "strong" . 

® This quantity is also called cross helicity (see [70]). 
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are intermittent not only in space, but also in time, the transport of turbulent 
energy described in CLV02b may be feasible. 

What is the density structure that we expect to see? 
First of all, we do not expect to see tight correlation between density and mag- 
netic field. Such sort of correlation is expected in the traditional static picture 
of the ISM. Introduction of turbulence in the picture of ISM complicates the 
analysis (see discussion in [99]; CLV02a). Our results confirm earlier claims (e.g. 
CLV02a; [85]) that magnetic field - density correlations may be weak. First of all, 
some magnetic field fluctuations arc related to Alfvenic turbulence which does 
not compress the medium. Second, slow modes in low /3 plasmas are essentially 
density perturbations that propagate along magnetic field and which marginally 
perturb magnetic fields. 

On the small scales we expect to see structures that anti-correlate with mag- 
netic field and caused by the new regime of turbulence below the ambipolar 
damping scale. We mentioned above that the simulations in Cho, Lazarian & 
Vishniac (CLV02c) and theoretical calculations in Lazarian, Vishniac & Cho [55] 
show that the magnetic field in a newly discovered regime of MHD turbulence 
can produce a shallow spectrum E(k) ^ kr^ spectrum of magnetic fluctuations. 
Calculations in CL03 suggest that this will translate in the corresponding shallow 
spectrum of density. For cold neutral medium (CNM; see Draine & Lazarian [18] 
for a list of idealized phases) the spectrum of density fluctuations can protrude 
from a fraction of parsec to a scale of ~ 100 AU (see CL03) . 

9.3 Astrophysical significance of fast modes 

Results of CL02 and CL03 show that the Alfven and slow modes exhibit Goldreich- 
Sridhar scale-dependent anisotropy. However, it would be very wrong to forget 
that fast modes are isotropic. It is possible to show that in many instances that 
difference makes fast modes very important. Consider two examples. 

Cosmic Ray Propagation. 
The propagation of cosmic rays is mainly determined by their interactions with 
electromagnetic fluctuations in the interstellar medium. The resonant interaction 
of cosmic ray particles with MHD turbulence has been repeatedly suggested as 
the main mechanism for scattering and isotropizing cosmic rays. In these studies, 
it is usually assumed that the turbulence is isotropic with a Kolmogorov spec- 
trum (e.g. [91]). Yan & Lazarian [101] identified fast modes as being responsible 
for cosmic ray scattering. The scattering by Alfvenic turbulence, which is the 
default for most of the theoretical constructions, is from 15 to 5 orders smaller 
than it is usually obtained using Kolmogorov model of Alfvenic turbulence (see 
Fig. 11). 

Dust Grain Dynamics. 
Turbulence induces relative dust grain motions and leads to grain-grain colli- 
sions. These collisions determine grain size distribution, which affects most dust 
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Fig. 11. Rate of cosmic ray scattering by Alfvenic {left panel) and fast (right Panel) 
modes of MHD turbulence. The scattering by Alfvenic modes is negligible, although it 
is substantially larger than an earlier estimate in Chandran [5]. The default for many 
researchers is Alfvenic turbulence with Kolmogorov spectrum (upper curve on the left 
panel). Fast modes scatter cosmic rays much more efficiently in spite of being partially 
damped in the ISM the dependence of scattering on plasma /3 is a new prediction. 
Prom Yan & Lazarian [101]. 



properties, including starlight absorption and H2 formation. Unfortunately, as 
in the case of cosmic rays, earlier work appealed to hydrodynamic turbulence 
to predict grain relative velocities. Lazarian &: Yan [56] and Yan & Lazarian 
[102] considered motions of charged grains in MHD turbulence and identified 
the direct interaction of the charged grains with fast modes as the principal 
mechanism for acceleration of grains with radius larger than 10~^ cm. Those 
modes can acceleration provide grains with supersonic velocities (see fig. 12). 

9.4 Other examples: from HII regions to gamma ray busts 

Lithwick & Goldreich [61] addressed the issue of the origin of density fluctuations 
within HII regions. There the gas pressure is larger than the magnetic pressure 
(the 'high /3' regime) and they conjectured that fast waves, which arc essentially 
sound waves, would be decoupled from the rest of the cascade. They found that 
density fluctuations are due to the slow mode and the entropy mode, which 
arc passively mixed by shear Alfvcn waves and follow a Kolmogorov spectrum. 
Our results in CL03 suggest that fast modes may be also an important source 
of density fluctuations. In addition, results on the new regime of turbulence 
(CLV02c, CL03) indicate that the new regime of turbulence can fluctuations on 
very small scales and this entails resumption of the turbulent cascade [55] that 
was not considered in [61]. 

The whole machinery of MHD turbulence scalings is required to deal with 
turbulence in gamma ray bursts (see [48]). There both fast and Alfven waves 
can transfer their energy to emitting electrons. However, the ways that they 
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Fig. 12. A new mechanism of grain acceleration is based on the direct interaction 
of charged grains with MHD turbulence. Gyroresonance of charged grains with fast 
modes dominates the acceleration of charged grains in Cold Neutral Medium (CNM). 
Within dark clouds the interaction of MHD turbulence and charged grains is important 
only for small grains. When grains decouple from magnetic field, Kolmogorov theory 
provides an OK estimate for grain velocities. Prom Yan & Lazarian [102] and Lazarian 
& Yan [56]. 

transfer their energy are different and this may result in important observational 
consequences. 

Heating of ISM and of Diffuse Ionized Gas (DIG), in particular, is another is- 
sue where imbalanced MHD turbulence is important (see CLV02a). Compression 
of molecular clouds by MHD turbulence (see [81]), stochastic magnetic rccon- 
nection [54,55] are other examples when it is essential to know the fundamental 
properties of compressible MHD. 

10 Summary 

In the paper, we have studied generation of compressible MHD modes. Wc have 
presented the statistics of compressible MHD turbulence for high, intermediate, 
and low /3 plasmas and for different sonic and Alfven Mach numbers. For sub- 
Alfvenic turbulence we provided the decomposition of turbulence into Alfven, 
slow and fast modes. We have found that the generation of compressible modes 
by Alfvenic modes is suppressed and, contrary to the common belief, the drain 
of energy from Alfven to compressible modes is marginal along the cascade. As 
the result the Alfvenic modes form a separate cascade with the properties sim- 
ilar to those of Goldreich-Sridhar cascade in incompressible media. As Alfven 
modes shear slow modes they impose their scaling on them. On the contrary, 
fast modes show isotropy for both magnetic- and gas-pressure dominated plas- 
mas. The new insight into compressible MHD entails important astrophysical 
consequences that range from the dynamics of star formation to the dynamics 
of gamma ray busts. 
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Appendix 

Let us consider a small perturbation in the presence of a strong mean mag- 
netic field. We write density, velocity, pressure, and magnetic field as the sum 
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of constant and fluctuating parts: p — > po + p, v — > vq + v, P — > Pq + and 

B Bq + b, respectively. We assume that vo = and that perturbation is 
small : p << po, etc. Ignoring the second and higher order contributions, we can 
rewrite the MHD equations as follows: 

|j+PoV-v = 0, (46) 

po-^+ V(aV) - —(V X b) X Bo = 0, (47) 

— + V X [v X Bo] = 0, (48) 

where we assume a polytropic equation of state: p = a? p with a? = jpo/po. We 
follow arguments in Thompson (1962) to derive magnetosonic waves. Let £,{r,t) 
be the displacement vector, so that d^/dt = v. Assuming that the displacements 
vanish at t = 0, we can integrate the equations as follows 

P + PoV • e = 0, (49) 
I = a2V(V • + (V X b) X Bo/47rpo, (50) 
b = V X (^ X Bo). (51) 

The momentum equation (eq. 50) becomes 

I = a2v(V • + [V X (V X (^ X Bo))] x Bo/47rpo 

= a2v(v . ^) + v(s2y . ^ _ Bo • VBo • 

- (Bo • Vf^/Anpo + [Bo (Bo • V)V • C] /47rpo (52) 
Using a = o? IV\ = l3{-f/2), Va = Bo/Anpo, we have 

Vi - V[(a + 1) V • ^ - (Bo • V) (Bo • 0] - (Bo • V)^^ + (Bo • V) (V • ^Bo = (53) 

In Fourier spacethc equation becomes 

l/yj + fck[(Q + l)kCk - + k^^ - kiikCkk = 0, (54) 

where Cfe = ^ " k, Cl| = C ' k = k/A;, and ky is unit vector parallel to Bo (i.e. 
k|| = Bo). Assuming ^ = —co^^ = —c^k'^^, we can rewrite (54) as 

{c^/Vi - cos^ 9)^ - [{a + l)a - cos^C||]k + cos^^^H = 0, (55) 

where cos6' = k\\/k and 6 is the angle between k and Bo- 
Using k = sin ^kj_ + cos ^k|| , we get 

(c^/Va - cos^ 0)^ - + 1)6 - cos 6*^11] sin 6'k_L 

-{[{a + l)^k - cos9^l\]cos0 - cose^k}^ = 0. (56) 
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Writing ^ = ^_Lk_L + ^||k|| + ^<^(p, we get 

(cVVa - cos^ e)^j_ - [{a + l)^k - cos^^ij] sin 61 = 0, 
(cV^i - cos^ 61)^11 - [a^k - cos6lC||] cos 61 = 0, 

(cV^i - cos2 e)^^ = 0. 



(57) 
(58) 
(59) 



The non-trivial solution of equation (59) is the Alfven wave, whose dispersion 
relation is uj/k ~ Vacos6'. The direction of the displacement vector for Alfven 
wave is parallel to the azimuthal basis (p: 



= -if = k± X k||. 



(60) 



Let us consider solutions of equations (57) and (58). Using = ^J-Sin^ + 

^11 COS0, we get 



(c^/Va - cos^ 0)U - (a + 1) sin^ 9^ - a cos 6 sin ff^n = 0, 
(cVv| - cos^ 6')^|| - a sin 6* cos 6*^^ - (a - 1) cos^ 6*^11 = 0. 
Rearranging these, we get 

(cV^i - a sin^ - 1)^_l - a cos 6* sin 6*^11 = 0, 
(cVVa - acos2 6')^|| - asin6'cos6'^_L = 0. 

Combining these two, we get 

(cVVj - a sin^ - \){^IVl - a cos" 6) 
= o? sin^ ^cos^ Q. 

Therefore, the dispersion relation is 

c^jVX - (1 + ay^lVl + a cos^ ^ = 0. 

The roots of the equation are 

1, 



4. = 



(1 + a) ± V(l + a?-4acos2^ 



where subscripts 'f and 's' stand for 'fast' and 'slow' waves, respectively. 
We can write 



/Ciikii + fc_Lk_L. 



^ = C||k|| +$_Lk_L oc 

Plugging eq. (67) into eq. (63) and (64), we get 

= a COS ^ sin 0^11 , 



1 + a . 2„ 1 

± a sm 0—1 

2 2 



1 + a , /D 2 

± a COS 

2 2 



^11 = acos0sin0^j_. 



(61) 
(62) 

(63) 
(64) 



(65) 
(66) 

(67) 
(68) 

(69) 
(70) 
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where D = {1 + a)^ — 4acos^ 6. Using k» = kcos6 and k± = kcos9, we get 



-1 



2 

1 + a 



± 



± 



2 2 
Arranging these, we get 



C-L^ll ~ asin^ ^^_lA;|| = acos^ ^^||fc_L, 
^\\k± — acos^ ^^ii^x = asin^ d^±k\\. 

^\\k± _-l + a±VD 



(71) 
(72) 

(73) 



where the upper signs are for fast mode and the lower signs for slow mode. 
Therefore, we get 



t(x{-l + a- /D)A;||k|| + (1 + a - VD)k±'k±, 
1/ cx (-1 + a + \/D)k\\k\\ + (1 + a + \/D)kA_kj_. 



(74) 
(75) 



The slow basis lies between k|| and —6. The slow basis ^/ lies between kj^ 
and k (Fig. 3). Here overall sign of and C/ is not important. 
When a ^ 0, equations (75) and (74) becomes 

~ kj| — (Q!sin6'cos^^)kJ^, (76) 
if « (asin6'cos6')k|| + k^. (77) 

In this limit, is mostly proportional to k|| and to kj_. When a — > oo, 
equations (75) and (74) becomes 

is ^ -9+ (sin e cos 6i/a)k, (78) 

^7 w (sin e cose/a)e + k. (79) 

When a = oo, slow modes are called psettdo- Alfvenic modes. 

We can obtain slow and fast velocity component by projecting Fourier veloc- 
ity component v^ onto and ^y, respectively. 

To separate slow and fast magnetic modes, we assume the linearized continu- 
ity equation {ujpk = pok • v^) and the induction equation (wbfe = k x (Bq x Vfe)) 
arc statistically true. From these, we get Fourier components of density and 
non-Alfvenic magnetic field: 



Pk = {poAvk,s/Cs)yi ■ is + {PO^VkJ / Cf)\<i ■ if 

= Pk,s + PkJ, 

hk = {BoAvk,s/cs)\^t> X is\ + (Bo^Wfc,//c/)|Bo x 
= bk,s + bkj, 

= Pfe,«(Bo/po)(|Bo X il/k-i) 
+ Pfc,/Wpo)(|Box|j|/k-|j), 



(80) 
(81) 
(82) 
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where Avk oc v'^ — v^^ (superscripts '+' and '-' represent opposite directions of 

wave propagation) and subscripts 's' and 'f stand for 'slow' and 'fast' modes, 
respectively. From equations (80), (81), and (82), we can obtain p^^s, Pkj, bk,sj 
and bkj in Fourier space. 



